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We report on our studies of fractal butterflies in biased bilayer graphene in the fractional quantum 
Hall effect (FQHE) regime. We have considered the case when the external periodic potential is 
present in one layer and have illustrated the effect of varying both the periodic potential strength 
and the bias voltage on the FQHE and the butterfly energy gaps. Interestingly, the butterfly 
spectra exhibits remarkable phase transitions between the FQHE gap and the butterfly gap for 
chiral electrons in bilayer graphene, by varying either the periodic potential strength or the bias 
voltage. We also find that, in addition to those phase transitions, by varying the bias voltage one 
can essentially control the periodic potential strength experienced by the electrons. 


Graphene §] placed on a hexagonal substrate with 
a twist is expected to display the Hofstadter butterfly 
pattern 0,0 in the energy spectrum of non-interacting 
Dirac fermions, when subjected to a perpendicular mag¬ 
netic field. The unique butterfly pattern was indeed dis¬ 
covered in 2013 in monolayer and bilayer graphene 0- 
0, placed on a hexagonal boron nitride substrate with 
rotational misalignment between graphene and the sub¬ 
strate. That arrangement resulted in the Moire pattern 
which actually introduces a large-scale periodicity in the 
Hamiltonian of the system, and the fractal butterfly pat¬ 
tern was the result of splitting of the Moire minibands 
(secondary Dirac cones) by the magnetic field that are 
exhibited in the magnetoconductance probe 0. Interest¬ 
ingly, a mathematical proof of the presence of the frac¬ 
tal pattern in the butterfly spectrum (the ‘Ten Martini 
Problem’) is also available in the literature 0- After 
that exciting experimental discovery of the fractal but¬ 
terflies, more recent studies (both theoretical [11] and 
experimental 0 ) have focused on the influence of the 
electron-electron interaction on the butterfly spectrum. 
Electronic properties of Dirac fermions in monolayer and 
bilayer graphene have been exhaustively studied in re¬ 
cent years [100. In a strong perpendicular magnetic 
field, interacting Dirac fermions [14j display the frac¬ 
tional quantum Hall Effect (FQHE) states [15] in mono- 
layer [lj| and bilayer graphene 0,0, that has also been 
experimentally observed [18]. The interaction effects in 
the fractal butterflies are, of course, more complex in the 
fractional quantum Hall effect regime, where one observes 
an interplay between the quantum Hall effect gap and the 
Hofstadter gap 0 |. Interestingly, in this work we find 
that the butterfly spectra exhibit remarkable phase tran¬ 
sitions for chiral electrons in bilayer graphene, where in 
addition to the phase transitions between the Hofstadter 
gap and the FQHE gap, one can essentially control the 
periodic potential strength experienced by the electrons 
by varying the bias voltage. 

We consider bilayer graphene with Bernal (AB) stack¬ 
ing in an external periodic potential with square symme¬ 
try 0, 0 - 0 . We label the layers of bilayer graphene 
by the indices 1,2 and assume that the periodic potential 


is present only in layer 1. Considering that the layers are 
stacked such that the B 1 and the A 2 sites are vertically 
aligned, the single-particle Hamiltonian of this system in 
a magnetic field (without the periodic potential) is writ¬ 
ten [1, 0 0, 0 
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two-dimensional electron momentum, A = (0, Bx, 0) is 
the vector potential, v F ~ 10 6 m/s is the Fermi ve¬ 
locity in graphene, U is the inter-layer bias voltage, 
7 1 ~ 0.4 eV is the interlayer hopping integral and 
£ = 1 for K valley and £ = —1 for K' valley. The 
corresponding wave function is described by a four- 
component spinor (i/j^ , i/j b , , "0b 2 ) T f° r valley K and 

0 for valle y K> 1 where and i’B 

are wave functions of sublattices A and B, respectively. 
We consider the fully spin-polarized electron system and 
therefore disregard the Zeeman energy. The eigenfunc¬ 
tion of the Hamiltonian m then has the form 
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where C 1 ,C 2 ,C 3 ,C i are constants and ip n j is the elec¬ 
tron wave function in the n-th Landau level (LL) with 
the parabolic dispersion, taking into account the periodic 
boundary conditions (PBC) 0,0. In the wave func¬ 
tion m the LL index n can take the values — 1,0,1,... 
and we assume that if the LL index of ip n 7 is negative 
then it is identically equal to zero. In this case, for 
n = — 1 the wave function © is T_ 1 ■ = (0, 0,0, ip 0 •) 
and there is only one energy level corresponding to this 
case. For = 0, C x = 0 and there are three energy states 
in this case. Following the convention for the indexing 
the energy levels introduced in [13], we label the states 
for n = — 1 and n = 0 as where i = —2, —1,1, 2 is 
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FIG. 1: The Landau levels for n = — 1 and n = 0 as a func¬ 
tion of the bias voltage U for two values of a, with periodic 
potential strength Vo = 20 meV. The numbers next to the 
curves denote the corresponding Landau level as described 
in the text. The regions where the gap corresponds to the 
FQHE (butterfly) gap are drawn as blue (red) dots, (a) and 
(c) correspond to the system with N e = 4 and (b) and (d) to 
N e = 6. 

the label of states in ascending order of the energy val¬ 
ues. In particular, for valley K the state corresponding 
to n = — 1 has the index 0^ and 0^ 1 for valley K'. 

The complete many-body Hamiltonian for this system 
can be written as 

N e N e 

n = Y^[nf + v (x i ,y i )\+-Y / v i j (3) 

i 

where the second term is the periodic potential, which is 
nonzero only for the components of layer 1 and the last 
term is the Coulomb interaction. We disregard the val¬ 
ley mixing terms (short-range interaction energies) due 
to the_ periodic potential and the Coulomb interaction 
[iiSl l25l |. Both the inter-layer bias voltage and the pe¬ 
riodic potential break spatial inversion symmetry and 
therefore the valley degeneracy is lifted in this system 
f2(|. For the many-body problem we consider a system 
of finite number N e of electrons in a toroidal geometry, 
i.e., the size of the system is L x = M x a 0 and L y = M y a 0 
( M x and M y are integers, a 0 is the period of the exter¬ 
nal potential) and apply periodic boundary conditions 
(PBC) in order to eliminate the boundary effects. Defin¬ 
ing the parameter a = 4> 0 /4> (the inverse of the magnetic 
flux through the unit cell measured in units of the flux 
quantum), where (f> = Ba q is the magnetic flux through 
the unit cell of the periodic potential and <f> 0 = he/e the 
flux quantum, we have N s /(M x M y ) = 1/a, where N s 
is the number of magnetic flux quanta passing through 
the system or, alternatively, it describes the LL degen¬ 
eracy for each value of the spin and valley index and 


FIG. 2: The dependence of both the FQHE and the butterfly 
gaps for N e = 4 electrons on the bias voltage U and the Lan¬ 
dau levels and 0j + ^ for two values of a, with the periodic 
potential strength Vo = 20 meV. The line corresponding to 
the FQHE (butterfly) case is depicted in red (blue). The re¬ 
gion of the bias voltage U where the gap corresponds to the 
FQHE gap is marked in cyan. 


takes integer values. The filling factor is defined as 
v = p/q = N e /N s , where p and q are again coprime in¬ 
tegers. In order to solve this problem we first construct 
the Hamiltonian matrix using the Hamiltonian operator 
© and the many-body states ■ ■ ■ iJn ) (besides j ?: , 

each single-particle state is characterized by the LL and 
valley indices which are not shown, but are implicitly 
assumed to be included in the indices j/) which are con¬ 
structed from the single-particle eigenvectors ©. After 
that we use the center of mass (CM) translation algebra 
PS 0, m S3 Hi and the eigenstates of CM transla¬ 
tions to bring the Hamiltonian matrix into block diag¬ 
onal form, where each block can then be diagonalized 
using the exact diagonalization procedure jl5j . 

In this work we consider the filling factor v = 1/3 
for a fully spin polarized system and for a = 1/2 and 
a = 1/3. The period of the external potential is taken to 
be a 0 = 20 ran, and the inter layer hopping integral to 
be 7 1 = 30 meV P3 , which can be achieved by applying 
an in-plane magnetic field pl|. In Fig. |T] the dependence 
of the single-particle LL 0- +) on the bias-voltage is pre¬ 
sented for the K valley. The color of the filled dots in¬ 
dicate the regions of the bias voltage where for N e = 4 
(a,c) and N e = 6 (b,d) the gap corresponds to the but¬ 
terfly region or the FQHE region. The periodic potential 
strength V 0 is taken to be V 0 = 20 meV. In order to 
understand the phase transitions observed in Fig. |T] the 
wave functions and also the impact of the bias voltage 
on these wave functions should be analyzed for each LL. 
As mentioned above, for level o(// the wave function has 
a non-zero component only in layer 2 and this remains 
true for all values of the bias voltage. Due to the fact 
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that the periodic potential is present only in the first 
layer, the wave functions and therefore also the FQHE 
gaps do not depend on the bias voltage for LL Oly and 
there is no phase transition in this case. In the LL 0^t,\ 
for U = 0 the electrons are mostly located in the layer 
2, although they have small probability of being in layer 
1. Increasing the bias voltage, both the single-particle 
and the many-particle system become even more polar¬ 
ized in layer 2, and therefore the periodic potential has a 
negligible impact on this level as well and we observe the 
FQHE gap for all values of U. The situation is different 
for LL 0^ and O^. For [7 = 0 the electrons in 
are mostly localized in layer 1 and therefore the periodic 
potential has a drastic impact in this case. 

In monolayer graphene the magnitude of the periodic 
potential V 0 = 20 meV for N e = 4 and N e = 6 and 
for both values of a, closeis the FQHE gap and opens 
the butterfly gap [ljj. It should be pointed out that for 
a = 1/2, there is no gap in the butterfly spectrum for 
the non-interacting electrons, but the electron-electron 
interaction opens a gap [Tlj| . However, for a = 1/3 the 
gap is due to both the Hofstadter gap observed in the 
single-particle case and due to the contribution from the 
electron-electron interaction. Therefore, the gap for LL 
0 A correspond to the butterfly gap for low values of 

the bias voltage U. The electrons in are mostly 
localized in layer 2 and therefore the periodic potential 
has only a minor effect on them. Hence for the 
LL, the FQHE gap is observed for low values of the bias 
voltage U. By increasing the bias voltage U there is 
an anticrossing between these two LLs (0^ and O^) 
and thereafter the layer polarization in each LL changes 
drastically. In particular, for the LL at the bias 
voltage U = 250 meV the probability of electrons being 
localized in layer 2 and for the electrons in LL to be 
localized in layer 1 are already « 0.95. This results in a 
phase transition in both LLs, namely the gap in the LL 
oJ + ' ) which initially represented the butterfly gap now 
correspond to the FQHE gap. The opposite behavior 
occurs for the LL O^. 

The closure of the FQHE gap by the external peri¬ 
odic potential also occurs in monolayer graphene 0- 
However, in bilayer graphene one can control the actual 
strength of the periodic potential experienced by the elec¬ 
trons essentially by applying the bias voltage. The impli¬ 
cations of this interesting result will be discussed below. 
The described behavior is almost the same for both values 
of a = 1/2 and a = 1/3. The essential difference between 
these two cases is that the butterfly gap for a = 1/3 is 
substantially bigger than that for a = 1/2, and there¬ 
fore the value of the bias voltage U, which is needed to 
observe the phase transition from the butterfly region to 
the FQHE region for the LL is bigger than that for 
a = 1/2. We have done similar studies for the K' valley 
and similar phase transitions were observed in that case 
as well for the LLs 0^ and in the negative energy 
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FIG. 3: The dependence of both the FQHE and the butterfly 
gaps for N e = 4 electrons on the periodic potential strength V 0 
for Landau levels 0^ and and for the bias voltages U = 0 
meV and U = 200 meV, and a = 1/2. The line corresponding 
to the FQHE (butterfly) is depicted in red (blue) color. The 
region of the periodic potential strength V 0 where the gap 
corresponds to FQHE gap is marked in cyan. 


region. 

In Fig. [2j the dependence of both the FQHE and the 
butterfly gaps for N e = 4 electrons on the bias voltage U 
for Landau levels 0^ and and for two values of a 
is shown for V 0 = 20 meV. The region of the bias volt¬ 
age U where the gap corresponds to the FQHE and the 
butterfly gap is also indicated. As was already pointed 
out, the layer polarization of the electrons changes dras¬ 
tically near the anticrossing point of the LLs and 
and the consequence of that can be clearly seen in 
the dependence of the gaps on the bias voltage U. In 
Fig. [2j in the regions further away from the anticrossing 
point the gaps are almost constant and fall rapidly to zero 
when approaching the anticrossing point. Also it can be 
clearly seen that the FQHE gap is almost the same for 
both LLs and for both values of a, whereas the butterfly 
gap is almost twice as big for a = 1/3 compared to that 
of a = 1/2 as explained above. 

The dependence of both the FQHE and the butter¬ 
fly gaps for N e = 4 electrons on the periodic potential 
strength V 0 for Landau levels and is shown in 
Fig.[3]for the bias voltages U = 0 meV and U = 200 meV, 
and for a = 1/2. The same dependence for a = 1/3 
is shown in Fig. [U The region of the periodic poten¬ 
tial strength V 0 where the gap corresponds to the FQHE 
and the butterfly gaps is also indicated. In Fig. [3] and 
IU similar phase transitions between the FQHE and the 
butterfly gaps is observed as well, although the depen¬ 
dence of the gap on the periodic potential strength V 0 
is almost linear in comparison to the dependence on the 
bias voltage U. No phase transition is observed for the 
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FIG. 4: Same as in Fig. [3]but for a = 1/3. 


LL and U = 0 meV, because as was noted above, in 
this case the electrons are mostly localized in layer 2 and 
the periodic potential has almost no impact on the phys¬ 
ical system. This feature is not observed for the case of 
U = 200 meV, because as we noted above, application of 
a bias voltage gradually polarizes the electrons for the LL 
from layer 2 to layer 1 and the effect of the periodic 
potential is apparent already for U = 200 meV. A similar 
behavior is observed for LL 0^ + \ where the application 
of the bias voltage results in widening of the FQHE gap 
region for both cases of a. While the butterfly gap is 
almost linear for a = 1/3 in Fig. HI which indicates that 
the main contribution here is due to the Hofstadter gap 
(single particle), for a = 1/2 (Fig. [3]) it deviates from the 
linear behavior in some cases. This indicates that the 
butterfly gap due to the electron-electron interaction is 
highly non-trivial. 

While the closure of the FQHE gap and opening of 


the butterfly gap due to the external periodic potential 
occurs in monolayer graphene fisij , one cannot control the 
periodic potential strength in that system. Therefore, 
there is no direct method to vary the periodic potential 
strength in the experiment. In bilayer graphene, as our 
present work indicates, this can be achieved by applying 
the bias voltage on the sample. Based on the observations 
above, variation of the bias voltage offers us the ability to 
control the polarization of the electrons between the two 
layers, which essentially translates to the control of the 
strength of the periodic potential, and at the same time 
exhibits the phase transition between the FQHE gap and 
the butterfly gap. This can have significant implications 
for experimental realization of the fractal butterflies in 
the FQHE regime. 

In conclusion, we have utilized the exact diagonaliza- 
tion scheme to study the FQHE and the butterfly gaps in 
bilayer graphene under the applied interlayer bias voltage 
and for the filling factor v = 1/3. We have considered 
the case when the external periodic potential is present 
in one layer and have illustrated the effect of varying 
both the periodic potential strength and the bias volt¬ 
age on the FQHE and the butterfly gaps. Two values of 
the parameter a were considered, namely a = 1/2 and 
a = 1/3. We found that by varying either the periodic 
potential strength or the bias voltage for some Landau 
levels in both valleys, a phase transition from the FQHE 
gap to the butterfly gap or vice versa can be observed. 
While the periodic potential strength is a characteristic 
of the sample used in the experiment and cannot be var¬ 
ied directly, our finding shows that by varying the bias 
voltage, change of the periodic potential strength actu¬ 
ally experienced by the electrons can be achieved, which 
can have a huge impact on the experimental investigation 
of the fractal butterflies in the FQHE region. 
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